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Abstract 

We propose a g-deformation of the sit(2)-invariant Schrodinger equa- 
tion of a spinless particle in a central potential, which allows us not only 
to determine a deformed spectrum and the corresponding eigenstates, as 
in other approaches, but also to calculate the expectation values of some 
physically-relevant operators. Here we consider the case of the isotropic 
harmonic oscillator and of the quadrupole operator governing its interac- 
tion with an external field. We obtain the spectrum and wave functions 
both for q £ R + and generic q 6 S 1 , and study the effects of the q- value 
range and of the arbitrariness in the su q (2) Casimir operator choice. We 
then show that the quadrupole operator in I — states provides a good 
measure of the deformation influence on the wave functions and on the 
Hilbert space spanned by them. 

Mathematics Subject Classifications : 81R50, 17B37. 
Keywords: Quantum Groups, Schrodinger Equation. 

1 Introduction 

Since the advent of quantum groups and quantum algebras |t|, |2]j| , there 
has been a lot of interest in deformations of the harmonic oscillator, since the 
latter plays a central role in the investigation of many physical systems. Most 
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studies were concerned with the one-dimensional oscillator j2[ ||, |T^|. Various 
g-deformed versions of standard quantum mechanics in the Schrodinger repre- 
sentation were proposed for the latter by using either the ordinary differentiation 
operator (see e.g. [jiTj), or a (/-differentiation one (see e.g. p6|, |l8|, |22|). 

The many-body problem for one-dimensional g-oscillators or, equivalently, 
the deformation of the TV-dimensional harmonic oscillator in cartesian coordi- 



nates was also considered f 1 , but relatively few works dealt with the deforma- 
tion of the same in radial and angular coordinates. In one approach |]l4j| , only 
the radial problem was deformed via a purely algebraic procedure based upon a 
so 9 (2, 1) so(iV) dynamical symmetry algebra. Only the N = 2 case was dealt 
with in detail, although extensions to higher dimension would in principle be 
feasible. 

In another type of analysis, use was made of a differential calculus on the 
./V-dimensional non-commutative Euclidean space Q to construct and solve a 
g-deformed Schrodinger equation for the isotropic harmonic oscillator [|[ [l0|, |], 
[HI . In such a setting, the underlying SO(A r ) symmetry gets replaced by an 
SO q (N) symmetry. In still another study the wave functions and energy 
spectrum of the g-isotropic oscillator were obtained in a rather indirect way 
from some previous results for the g-linear oscillator |^| . The latter were derived 
by replacing classical Poisson brackets by g-commutators (instead of standard 
commutators). In the last two approaches, the wave functions involve some 
non-commutative objects: the variables of the quantum Euclidean space in the 
former, and the elements of the SU g (2) quantum group in the latter. This 
complicates the calculation and interpretation of operator matrix elements. 

The purpose of this Letter is to present an entirely different approach to 
the g-deformation of the su(2)-invariant Schrodinger equation of a spinless par- 
ticle in a central potential, which allows us not only to determine a deformed 
spectrum as in other works, but also to easily calculate the expectation values 
of some physically-relevant operators, thereby evaluating the deformation influ- 
ence on the corresponding wave functions. Here we consider in detail the case 
of a particle in an isotropic harmonic oscillator and of the quadrupole operator 
governing its interaction with an external field. 

In our approach, only the angular sector is deformed by using a represen- 
tation of the su g (2) quantum algebra on the two-dimensional sphere ^0|. This 
gives rise to an appropriate change in the angular part of the scalar product |l3| , 
and to the substitution of the su q (2) Casimir operator eigenvalue for the su(2) 
one in the radial Schrodinger equation fl2]| . The latter step may be performed 
in various ways since there is no unique rule for constructing the su 9 (2) Casimir 
operator. Similarly, the deforming parameter q may be assumed either real and 
positive, or on the unit circle in the complex plane (but different from a root 
of unity), provided different scalar products are used We will study the 

effects of these two choices on the solutions of the radial Schrodinger equation. 

In Section^], the su g (2)-invariant Schrodinger equation of the three-dimensional 
harmonic oscillator is introduced and solved. In Section |3|, its spectrum is stud- 
ied in detail for various choices of su g (2) Casimir operators and q ranges. The 
effect of the deformation on the corresponding wave functions is determined 
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in Section |4| by calculating the quadrupole moment in I = states. Finally, 
Section contains the conclusion. 



2 sw g (2)-Invariant Schrodinger Equation 

Let 

h 2 ( d 2 2 d C q \ 1 2 2 

= - ( ^-o + r- - dr ) + o^ 2 (i) 



2/i \<9r 2 r <9r r 2 J 2 

be the Hamiltonian of a g-deformed three-dimensional harmonic oscillator in 
spherical coordinates r, 8, (f>. Here C g is the su q {2) Casimir operator, which we 
may take as 

,12 

7. ( 2 ) 



C 9 = J+J- 



*"5 



where [x] ? = (q x - q~ x ) / (g - g" 1 ), and g = e" 1 e R+ or g = e™ G S* 1 (but 
different from a root of unity). The operators J3, J + , J_, satisfying the su q {2) 
commutation relations 

[J 3 ,J±]=±J±, [J+, J_] = [2J ]„ (3) 

are defined in terms of the angular variables by 

J3 = -i<?0, 

J+ = -j* (tan(0/2)pi], 9 T2 + cot(0/2) g Tl [T a ],) , 

J- = e- 1 *(cot(fl/2)[Ti],« T »+taii(fl/2)« Tl [T 2 ],) J (4) 

with Ti = -| (sin(9<9 e - ify), T 2 = - i (sin(9<9 e + ify) §f| @. 
Instead of Equation (|^), we may alternatively use the operator 

C' q = J+J- + [.h] q [.h~l] q (5) 

in Equation (|l|), in which case the corresponding Hamiltonian will be denoted 
hjH' q . 

The Hamiltonians H q and H' q remain invariant under su q (2) since they com- 
mute with J3, J + , J_, and they coincide with the Hamiltonian of the standard 
three-dimensional isotropic oscillator when q = 1. For simplicity's sake, we shall 
henceforth adopt units wherein h = fJ,= lo = 1. 

The su g (2)-invariant Schrodinger equation 

-Hg 4>nlmq(r, 0, 4>) = E nlq lp n l m q{r, 0, (/)) (6) 

is separable and the corresponding wave functions can be written as 

4>nhnq (V, 9,<f>) = H n lq [r)Yi mq (8,<f>). (7) 
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The q-spherical harmonics Yi mq (9,(f>) 1 satisfying the equations 

1 1 2 

C q Y lmq (e,0) = C q (l)Y lmq (8,<t>), C q (l) = 



J3Ylmq(0,(t)) = mY"; m? (6>, 0), 

are given by [§0|, [u| 

Y lrnq (0, 4>) = Nim q Qi q (cot 2 (6/2)) R l m0q (cot 2 (9/2)) cot" 1 (0/2) e im *, 

M mq = (-1)' ([2/ + l] q [l + m] q !) 1/2 (4tt[Z - m} q \y 1/2 , 
Q lq (cot 2 (tf/2)) = ni= (1 + q 2k - 21 cot 2 (0/2))~ 1 , 

R\ Mq (cot 2 (fl/2)) = [l] g \ [Z - m] q \ £ fe CT^feS^W^T ' 



(8) 
(9) 



(10) 



and [a:],! = [x] q [x- l] q . . . [1], if a; € 7V+, [0]„! = 1, and ([ar] 5 !) 1 = if x e Af~ . 
The functions Yi rnq (6, 4>), where I = 0, 1, 2, ... , and m = — Z, — Z + 1, . . . , Z, 
form an orthonormal set with respect to the scalar product |13|] 



(l'm'\lm) q = (q-q- 1 )(4lnq)- 1 £ d6 sin<9 J ( 
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sin 2 (9/2)+g- 2 cos 2 (0/2) 



q sineds - l Y lmq (6,<t>) 



(ii) 

where the upper (resp. lower) signs correspond to q € R + (resp. q € 5 1 ). Note 
that the definitions (|j) and (10) slightly differ from those given in ^0|. With 
the present choice, in the q — *■ 1 limit they go over into the definitions of su(2) 
generators and spherical harmonics used in most quantum mechanics textbooks 
(see e.g. §). 

The radial wave functions 1Z n i q (r) = r 1 S n i q (r) in Equation (0) are the 
solutions of the radial equation 



^-^ 1 -r 2 + 2E nlq y nlq (r) = 
that satisfy the condition : 

-WO) = 0, 

and are square integrable with respect to the usual scalar product, i.e. 



(n'l\nl) q 



dr S n 'i q (r)S n i q (r) = S n ^ n . 



(12) 
(13) 

(14) 



They are given by 



S n i q (r) = v/2p) T [ai + n + 



-1/2 



e~2 r r ai L" 



? (r 2 ), (15) 
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where L„ l 2 (r 2 ) is an associated Laguerre polynomial pL and ai is any one of 
the two solutions 



a l± = l -±\ q {l), X q {l) = ^\+C q {ll (16) 

of the equation 

ai(ai - 1) = C,(0, (17) 

provided it satishes the condition ai S The corresponding energy eigenval- 
ues are 

E n i q = 2n + ai + -, n = 0,1,2,.... (18) 

The ^(N + 1)(N + 2)-degeneracy of the isotropic oscillator energy levels, where 
N = 2n + I, is therefore lifted. 

Similar results hold for the choice (||) for the Casimir operator, the only 
change being the substitution of C' q {l) = [l] q [l + l] q for C q (l). To distinguish 
the latter choice from the former, we shall denote all quantities referring to it 
by primed letters (X' q (l), a' v E' nlq , K lq (r), S' nlq {r), ...). 



3 Spectrum of the su g (2)-Invariant Harmonic 
Oscillator 

In the present section, we will study the condition ai E R + for the existence of 
the radial wave functions (|l5|) and of the corresponding energy eigenvalues (|l8|) , 
as well as the behaviour of the latter as functions of I and q for the two choices 
(||), ((H) of Casimir operators, and for q E R + or g e S 1 . 

Let us first consider the case where q = e w E R + . Since the spectrum is 
clearly invariant under the substitution q — > q^ 1 , we may assume q > 1, i.e., 
w > 0. 

In the C q case, A?(Z) = sinh 2 ((Z + \/2)w)/ sinh 2 w > for I = 0, 1, 2, . . . , 
hence both roots ((l^) of Equation (|l7]) are real and distinct. However, for 
/ 7^ 0, X q (l) = sinh((Z + l/2)u;)/ sinh w > 1, showing that only ai + is positive 
and therefore admissible, whereas for I — 0, A g (0) = [2 cosh(w/2)] _1 < 1/2 if 
q 7^ 1, so that ao+ and ao- are both admissible. Note that in the undeformed 
case (q = 1), one gets Ai(0) = 1/2, so that the root ao- has then to be discarded 
in accordance with known results. For q ^ 1, the spectrum therefore comprises 
the energy eigenvalues 

E n0q± = 2n + l±- ±— -, (19) 

2cosh(u)/2) 

sinh((/ + 1/2) W ) 

£;„i g = 2n+lH — , t = 1, 2, (20) 

sinn w 

The appearance of an additional I — level was already observed in the deformed 
Coulomb potential case p~2|. 
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In the C' q case, \' q 2 {l) = [l] q [l + l] g + \ > \ for I = 0, 1, 2, . . . , so that 
both roots a' l+ and a[_ are real and distinct, but only the former is positive, 
hence admissible. The spectrum therefore comprises the same levels as in the 
undeformed case, their energies being now 



[4sinh(^) sinh((Z + l)w) + sinh 2 w] 1 / 2 
2 sinh w 

Note that the energy of the / = states is left undeformed : 



^mi^ sillily -r ±jwj -r amii w - 

£„z 9 = 2n + l + o > i = 0,1,2,... (21) 



E' n0q = E n0 = 2n+~. (22) 

Expanding the right-hand sides of Eqs. (p9|), (|2C|), and (^TJ) into powers of 
w shows that in the neighbourhood of q — 1, i.e., w — 0, 

£„ 0g± ~2 II + 1±1t1 W 2 (1-| W 2 + ..-), 

£„ ig ~2n + Z + § + 3 L(2Z- 1)(2Z + l)(2i + 3)w 2 

x {1 + 3^[12l(i + 1) - 25] W 2 + •••}, I = 1, 2, ... , 

(23) 

+ ' + l + ilTiT- 2 {2^ + l)-l 

24i 3 (i+l) 3 -56i 2 (i+l) 2 -10Z(i+l)+7 2 , \ 7 _ n 1 9 

60[4Z(Z+1) + 1] ^ I' ' ' ' 

Hence, for I ^ 0, -E n ; 9 and are increasing functions of w in the neighbour- 
hood of w — 0, whereas for Z = 0, E n0q+ and E n0q _ have opposite behaviours, 
while E^Qg is independent of w. Moreover, for a given w value, the influence of 
the deformation increases with Z. Such trends are confirmed by Figure 1, where 
the first few lowest eigenvalues E$i q are plotted in terms of w. For the Z and w 
values considered, E' olq cannot be distinguished from E i q for Z ^ 0, or E o q + 
for / (I. 

Let us next consider the case where q — e lw G S 1 . Owing to the invariance of 
the spectrum under the substitution q — > we may now assume < w < n. 

For C q , A 2 (Z) = sin 2 ((Z + l/2)i(j)/sin 2 w vanishes for w = kn/(l + l/2), k = l, 
2, . . . , Z, but these w values are in any case excluded as roots of unity. Hence, 
for the w values considered here, Equation (|l7]) has two real, distinct roots. If 
\ q (l) = |sin((Z + l/2)w)/ smw\ < 1/2, both roots a; + and a/_ are admissible, 
whereas if X q (l) > 1/2, only ai + is so. The conditions can be reformulated in 
terms of 

7,(0 = 4sin 2 iz;C* g (Z) = -[-4cos((2Z + 1) w) + cos2w + 3]. (24) 

One finds two admissible roots if 7 g (Z) < 0, but only one root a/+ if 7 9 (Z) > 0. 
For instance, for Z = 0, there is a single eigenvalue for any w value, 

E n0q = 2,i + l + l —— , (25) 

2cos(iz;/2) 
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Figure 1: Spectrum of n = states in terms of w for q = e w e R + . The solid 
lines correspond to E 00q+ and E i q , I = 1, 2, . . . , 6, and the crossed line to 

EoOq- ■ 
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while for 1 = 1, there are either two or one eigenvalues, which are given by 



E nlq± =2n+l± 4 -Jg / y , if^lI<cos W <-^, 
Eni q -2n+l+ 4 -;gy , if cosw (26) 

or cosw > ~ 7 \^ , 

respectively. 

By proceeding in a similar way for C' q , one finds in terms of 

i q {l) = 4sin 2 wC' q (l) = 4sin((Z + l)w)sm(lw) (27) 

that there are two admissible roots a' l+ , a[_ if — sin 2 w < t£(Z) < 0, only one 
a' l+ if either "f' q {l) > or "f' q {l) — — sin 2 w (in which case a' l+ = 1/2), or none 
if j' q (l) < — sin 2 w. For instance, for I = 0, there is a single eigenvalue fl2^ ) for 
any w value and it coincides with the undeformed one , while for 1 = 1, there 
are two or one eigenvalues if cosw > —1/8, 

E' nl ± = 2n + 1 ± \\/l + 8 cosw, if — | < cosw < 0, 

(28) 

E' nlq = 2n + 1 + ~\/T+~8 cosw, if cosw = — g or cosw > 0, 

and none if cosw < —1/8. 

In both the C q and C' q cases, similar results can be derived for higher / 
values, and close enough to w = 0, for given n and I values one always finds 
a single eigenvalue going into the undeformed one, E n i = 2n + I + 3/2, for 
q — > 1. For all / values, the expansion of this eigenvalue into powers of w can be 



obtained from Equation (23) by substituting iw for w. Hence, in a small enough 
neighbourhood of w = 0, for / ^ 0, E n i q and E' nlq are decreasing functions of w 
, whereas E n o q is increasing and E' n0 remains constant. 

The first few eigenvalues Eoi q , going into Eqi when q — * 1, are displayed on 
Figure 2 for < w < it. One should remember that for some I and w values, 
there may exist other eigenvalues, which are not plotted on the figure, and that 
the discrete set of points, where w is equal to a root of unity, is excluded. Such 
is the case, in particular, of the points where Eoi q (I ^ 0) takes its minimal 
value 1, and of the point w = it/2, where Eoi q = 1 + (l/y/2) for any I value. 
The influence of the deformation on the spectrum is rather striking. For high w 
values, the levels get mixed in a very complicated way. It is remarkable that in 
the neighbourhood of w = tt/2, one obtains a spectrum very close to one with 
equidistant, infinitely-degenerate levels. 

In the C' q case, the situation is still more complex as some I values may 
disappear on some intervals. In the neighbourhood of w = 0, however, as in the 
real q case, E' nlq cannot be distinguished from E n \ q . 
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Figure 2: Spectrum of n = states in terms of w for q = e lw G 5 1 . The solid 
and dashed lines correspond to Eoi q (or Eoi q+ if there exists another eigenvalue 
with n — and the same I), where I = 0, 2, and Z = 1,3, respectively. 
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4 Quadrupole Moment in / = States 



The purpose of the present section is to study the effect of the deformation on the 
wave functions of the su g (2)-invariant harmonic oscillator, given in Equation ([?]), 
(|l0|), and ([To]), by determining the variation with w of the expectation value 
of some physically-relevant operator. For the latter, we choose the electric 
quadrupole moment operator, and we consider the quadrupole moment in a 
state with definite n, I values, which is defined conventionally as Q 

Q nlq = (nll\(3z 2 ~r 2 )\nll) q . (29) 

Equation ( p9[ ) corresponds to the choice C q for the su q (2) Casimir operator. 
When using instead C' q , the quadrupole moment will be denoted by Q' n i q - 
The undeformed counterpart of Q n i q and Q' n i q is given by 

Q nl = (nl\r 2 \nl) (3 cos 2 9 - l))\ll) = Un + l + ^\ f-^Tg) ■ ( 30 ) 

It vanishes for I = as a result of the familiar selection rule for the angular 
momenta I and 2 coupling. Nonvanishing values of Q„0q or Q' n0q will therefore be 
a direct measure of the effect of the deformation. Note that since in the C q case, 
there are two energy eigenvalues with I — and a given n value for q E R + (see 
Equation (|l9|)), we have to distinguish the corresponding quadrupole moments 
by a ± subscript. 

As in Equation (|30|), Q n oq can be factorized into radial and angular matrix 
elements, 

Qnoq = (n0|r 2 |n0) 9 (00| (3 cos 2 9 - l))|00),. (31) 

For any / value, the former is simply obtained by replacing I by a/ — 1 in the 
undeformed radial matrix element. Hence 

(n0\r 2 \n0) q = 2n + a + ^. (32) 

The calculation of the latter is more complicated as it implies the use of the 
deformed angular scalar product (|IT|), 

(OO|(3cos 2 - 1)|00), - 3( 1 g ~ 1 g " 1) (2, + V0 - I- ( 33 ) 

where 

1„= dOsmO d(/) 5 q sin6do-l 2 g / 34 x 

9 Jo Jo sin 2 (9/2)+ q- 2 cos 2 (9/2) H V ' 

This integral can be easily performed by making the changes of variables z = 
pe^, z = pe'^, p = cot(6»/2), and r) = p 2 @, |l|. One gets 

An [°° (l-q- 2 y) 2 q + q- 1 ( q + q~\ \ f . 
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Introducing Equation (pq) into Equation (|33|), we finally obtain 



Hence 

(00| (3 cos 2 0-1)100), = 2c g TnhC +1 - S. if? = e-Gi?+, 

2 cos 2 w-\-l _|_ 3 cos w if (7 F 5^ 



(37) 



The deformed quadrupole moment Q n oq is therefore an even function of w, so 
that we may again restrict ourselves to0<w<ooor0<w<7r according to 
whether q is real or complex. 



The deformed radial matrix element (32) (or its counterpart for Q' n Q„), being 



just equal to the corresponding energy eigenvalue in the units used, varies in 
the same way with w. From Equations (|l9|), ( |2^ ) and (^5|), it follows that for 
increasing w, it decreases (resp. increases) from 2n + | (resp. 2n + 5) to 2n + 1 
for C 9 , g e and «o+ (resp. Q!o_), increases from 2n + | to +00 for C q , 
q E S 1 , and remains constant for C' q , q E R + or q E S 1 . The effect of the 
deformation is therefore not significant, except in the case of C q when q E S 1 . 

For real q, the deformed angular matrix element, given in Equation (^7|), 
increases from to 2 when w goes from to +00. Hence, it is obvious that 
QnQq- and Q'nQq are increasing positive functions of w. It can be checked that 
the same is true for Q n oq+, the variation of the angular part of the matrix 
element dominating that of the radial one. On Figure 3, Q n oq+, QnOq-, and 
Q' n0q are displayed in terms of w for n = and n = 1. One can see some effect of 
the choice of Casimir operator and ao root, but it becomes significant only for a 
very large deformation. It should be stressed that the undeformed quadrupole 
moments in I ^ states and the deformed ones in Z = states have opposite 
signs. Both become comparable in absolute value for w > 1. 

For complex q, the deformed angular matrix element ( |37| ) decreases from 
to — 00 when w goes from to n. Hence, both Q n oq and Q' n o q are decreasing 
negative functions of w. On Figure 4, they are displayed in terms of w for n = 
and n = 1. Apart from the sign, which is now the same as that of Q n i for I 7^ 0, 
the conclusions remain similar to those for the real q case. 



5 Conclusion 

In the present Letter, we did show that as those of the free particle and of the 
Coulomb potential [[I2| , the su g (2)-invariant Schrodinger equation of the three- 
dimensional harmonic oscillator can be easily solved not only for q E R + , but 
also for generic q E S . It is worth stressing that we have been working in 
the framework of the usual Schrodinger equation (i.e., with no non-commuting 
objects contrary to some other approaches H, |l(| [|, |l9|, but with wave 
functions belonging to a Hilbert space different from the usual one, since the 
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Figure 3: Quadrupole moment of I = states in terms of w for q = e w e R + 
The solid, crossed, and dashed lines correspond to Q n oq+, QnOq-, and Q' n0q 
respectively, for n = (three lowest curves) and n = 1 (three highest ones). 
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Figure 4: Quadrupole moment of I = states in terms of w for q = e™ e S 1 . 
The solid and crossed lines correspond to Q n oq an d Q' n0q , respectively, for n = 
(two highest curves) and n = 1 (two lowest ones). 



13 



angular part of the scalar product has been modified when going from su(2) to 
su q {2) @. 

In the real q case, we did show that the spectrum is rather similar to the 
undeformed one, except that the energy levels are no more equidistant and that 
their degeneracy is lifted. For a given n value, the spacing between adjacent 
levels corresponding to I and I + 1, respectively, increases with / and with the 
deformation. In addition, there appears a supplementary series of I = levels 
when the Casimir operator C q is used. Apart from this, for small deformations, 
the results are rather insensitive to the choice made for the Casimir operator. 

In the complex q case, we did show that the spectrum is more complicated 
as for I 7^ and any n value, there may exist 0, 1, or 2 levels according to 
the deformation. The existence or inexistence of levels is also rather sensitive 
to the choice made for the Casimir operator. Close enough to q — 1, there 
however always exists a single level going into the undeformed one for q — > 1. 
In that region, the spacing between adjacent levels corresponding to I and l + l, 
respectively, now decreases with / and with the deformation. 

The closeness of our approach to the standard one did also allow us to 
study the effect of the deformation on the wave functions and the Hilbert space 
spanned by them. We did establish that it is rather strong as the quadrupole 
moment in the I = states, which vanishes in the undeformed case, now assumes 
a positive (resp. negative) value for q 6 R + (rcsp. q e S 1 ) irrespective of the 
Casimir operator used. 
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